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Abstract 

Using the influence functional formalism we show how to derive a generalized Ein- 
stein equation in the form of a Langevin equation for the description of the backreaction 
of quantum fields and their fluctuations on the dynamics of curved spacetimes. We 
show how a functional expansion on the influence functional gives the cumulants of the 
stochastic source, and how these cumulants enter in the equations of motion as noise 
sources. We derive an expression for the influence functional in terms of the Bogol- 
ubov coefficients governing the creation and annihilation operators of the Fock spaces 
at different times, thus relating it to the difference in particle creation in different his- 
tories. We then apply this to the case of a free quantum scalar field in a spatially flat 
Friedmann- Robertson- Walker universe and derive the Einstein-Langevin equations for 
the scale factor for these semiclassical cosmologies. This approach based on statisti- 
cal field theory extends the conventional theory of semiclassical gravity based on a 
semiclassical Einstein equation with a source given by the average value of the energy 
momentum tensor, thus making it possible to probe into the statistical properties of 
quantum fields like noise, fluctuations, entropy, decoherence and dissipation. Recog- 
nition of the stochastic nature of semiclassical gravity is an essential step towards the 
investigation of the behavior of fluctuations, instability and phase transition processes 
associated with the crossover to quantum gravity. 
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1 Introduction 



Backreaction of quantum processes like particle creation in cosmological spacetimes |l] has 
been considered by many researchers in the past for the purpose of understanding how 
quantum effects affect the structure and dynamics of the early universe near the Planck time 
@, |J. Because of the general nature and complexity of the problem, backreaction studies 
have also been used as a testing ground for the development and application of different 
formalisms in quantum field theory in curved spacetime Q, e.g., regularization schemes to 
obtain finite energy-momemtum tensor, perturbation methods, effective action formalism, 
etc. The most recent stage of development for the discussion of cosmological backreaction 
problems was the use of Schwinger-Keldysh (or closed-time-path, CTP) functional formalism 
||, which, being formulated in the in — in boundary condition, gives rise to a real and causal 
equation of motion ( the semiclassical Einstein equation), where the expectation value of 
the energy-momentum tensor of a quantum field acts as a source which drives the classical 
effective geometry. In this equation one can identify a nonlocal kernel in the dissipative term 
whose integrated dissipative power has been shown to be equal to the energy density of the 
total number of particles created, thus establishing the dissipative nature of the backreaction 
process || [7j . 

In pursuing a deeper understanding of the statistical mechanics meaning of dissipation, 
one of us |5J cast this backreaction problem into the conceptual framework of quantum open 
systems ||. He made the observation that a Langevin-type equation is what should be 
expected, and predicted that for quantum fields a colored noise source should appear in the 
driving term. He also conjectured that the particle creation backreaction problem can be 
understood succintly as the manifestation of a general fluctuation-dissipation relation for 
quantum fields in dynamical spacetimes, a non-equilibrium generalization of such relations 



depicting particle creation in black holes [yj], [11] and de Sitter universe |12|1 . The missing 
piece in this search is the noise term associated with quantum fields. 

To look into this aspect of the backreaction problem in semiclassical gravity, as well as 
exploring the quantum origin of noise and fluctuations in inflationary cosmology [TBI, and 



understanding the decoherence problem in quantum to classical transition jlJJ , Hu, Paz and 



Zhang II)] looked into the relation of colored noise and nonlocal dissipation in a quantum 



Brownian motion model with the influence functional of Feynman and Vernon [17], [19 
In this formalism the effects of noise and dissipation can be extracted from the noise and 
dissipation kernels in the real and imaginary parts of the influence functional, their interrela- 
tion residing in the fluctuation-dissipation relation obtained as a simple functional relation. 
If one views the quantum field as the environment and spacetime as the system in the quan- 
tum open system paradigm, then the statistical mechanical meaning of the backreaction 
problem in semiclassical cosmology can be understood more clearly M. In particular, one 
can identify noise with the coarse-grained quantum fields |2fj, |21j, derive the semiclassical 
Einstein equation as a Langevin equation pifl , and understand the backreaction process as 



the manifestation of a fluctuation-dissipation relation [23|. Continuing their investigation of 



the backreaction problem via the CTP formalism, Calzetta and Hu [2"2j also found that the 
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results obtained from the influence functional formalism is the same as that obtained earlier 
(but displayed only partially) from the Schwinger-Keldysh method. This paradigm has also 
been applied to problems in quantum cosmology [p^|. (For an account of the search and 
discovery of these ideas, see |^5| p6|.) 

The specific goal of this paper is to derive the semiclassical Einstein equation in the form 
of a Langevin equation. Our primary task is the derivation of noise from the quantum field 
source, and we do this by carrying out a cumulant expansion of the influence functional. This 
goal is shared by two other companion papers addressing different aspects of this problem: 



In Ref. HU, using the closed-time-path method || 0, Calzetta and Hu identified the source 
of decoherence and particle creation to the noise kernel and showed their relation through 
the Bogolubov coefficients. They also showed the relation of quantum noise with classical 



fluctuations, and derived the semiclassical Einstein equation with a noise term. In Ref. 
Hu and Sinha started with the density matrix of the universe in quantum cosmology in the 
manner of [2~I and demonstrated the existence of a fluctuation-dissipation relation for the 
particle creation and backreaction problem in a Bianchi Type-I universe. These two papers 
together with this one present a quantum open system approach to the backreaction problem 
in semiclassical gravity and cosmology. Together they can serve as a platform for exploring 
the transition to quantum cosmology. It can also address the dissipative nature of effective 
theories || p7f , and (to the extent that Einstein's general relativity can be viewed as an 
effective theory) possible dissipative effects in the low-energy limit of any theory of quantum 



gravity. For a general discussion of these ideas, see p8 



This paper is organized as follows: In Sec. 2 we give a brief review of the influence 
functional formalism, mainly to establish notations. Readers familiar with it can skip to the 
next section. In Sec 3 we show how a functional expansion on the influence functional gives 
the cumulants of the stochastic source, and how these cumulants enter into the equations 
of motion as noise sources. In Sec. 4, following Plj, we derive the form of the Hamiltonian 
for a scalar field in terms of its normal modes and consider a class of actions where the field 
modes are coupled parametrically to the scale factor of the universe. We derive an expression 
for the influence functional in terms of the Bogolubov coefficients governing the creation and 
annihilation operators of the Fock spaces at different times which signify particle creation. 
In Sec 5, we present two standard cases of cosmological particle creation and derive the 
Einstein-Langevin equations describing its backreaction on the background spacetime. 



2 Influence Functional Theory 

Consider the quantum system described by the action 

S[a,q] = S[a] + S e [q\ + S int [a,q). (2.1) 

We will consider a to be our system variable and q to be our environmental variables. 
Typically the environment has infinite degrees of freedom which is denoted here by a bold 
type. 
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We will briefly review here the Feynman- Vernon influence functional method for deriving 
the evolution operator. The method provides an easy way to obtain a functional representa- 
tion for the evolution operator J r for the reduced density matrix p r . Let us start first with 
the evolution operator J for the full density matrix p defined by 

p(t) = J(t,t i )p(t l ). (2.2) 

As p evolves unitarily under the action of (2.1), the evolution operator J has a simple 
path integral representation. In the position basis, the matrix elements of the evolution 
operator are given by 



J(a f ,q f ,a' f ,q! f ,t | a h q i: a-, q-,tj) = tC(a f ,q f ,t | a i: q i: ti))C*(a' f , q' f , t | a-,q-,^) 

a f 1/ . a 'f q '/ . (o o\ 

t r % f f i \ ) 

= J Da J Dq exp — S[a, q] J Da! J Dq exp — — S[a', q'} 

ai q,: a [ q( 

where the operator /C is the evolution operator for the wave functions. In a path-integral 
representation, the functional integrals are over all histories compatible with the boundary 
conditions. We have used the subscripts i, f to denote the initial and final variables at ti,t. 
The reduced density matrix is defined as 



^oo +oo 



p r (a,a')= J dq J dq'p(a, q; a', q')5(q - q') (2.4) 

— oo — oo 

and is propagated in time by the evolution operator J r 

+oo +oo 

p r (a,a',t) = / da,i I da[ J r {a,a! \t \ ai^a^ti) priai^al^ti ). (2-5) 



-oo — oo 



By using the functional representation of the full density matrix evolution operator given 
in (2.3), we can also represent J r in path integral form. In general, the expression is very 
complicated since the evolution operator J r depends on the initial state. If we assume that 
at a given time t = ti the system and the environment are uncorrelated 

pit = = p s {ti) x p e (ti), (2.6) 

then the evolution operator for the reduced density matrix does not depend on the initial 



state of the system and can be written [17 



as 



J~ r (af,a'f,t | ai^a'^ti) = J Da J Da exp— ^[a] — 5[a']| J 7 [a, a'} (2.7) 



4 



The factor jF[a, a'], called the 'influence functional', is defined as 



-oo +00 If If 



jF[a,a'] = J dq_f J dq { J dq[ J Dq J Dq 



—00 —00 



ir 1 (2-8) 

x exp - [S e [q] + S int [a, q] - S e [q'\ - 5 int [a', q'])p e (qi, q-, U) 

= exp ^-S'/^fa, a'] 

where Sjf[cl, a'} is the influence action. The effective action for the open quantum system is 
defined as S e ff[a, a'} = S[a] — S[a'} + SVfK a']- 

It is not difficult to show that (2.8) has the representation independent form [21| 



F[a,a'} = Tr(u[a t , u }p e (U)W[a' t ^}) (2.9) 

where U is the quantum propagator for the action S e [q] + S int [a(s), q] where a(s) is treated 
as a time dependent classical forcing term. We have found this form to be very convenient 
for deriving the influence functional. 

It is obvious from its definition that if the interaction term is zero, the influence functional 
is equal to unity and the influence action is zero. In general, the influence functional is a 
highly non-local object. Not only does it depend on the time history, but -and this is 
the more important property- it also irreducibly mixes the two sets of histories in the path 
integral of (2.7). Note that the histories a and a' could be interpreted as moving forward and 
backward in time respectively. Viewed in this way, one can see the similarity of the influence 



functional |17| and the generating functional in the closed-time-path (CTP or Schwinger- 
Keldysh) integral formalism 0. The Feynman rules derived in the CTP method are very 
useful for computing the IF. 

In those cases where the initial decoupling condition (2.6) is satisfied, the influence func- 
tional depends only on the initial state of the environment. The influence functional method 
can be extended to more general conditions, such as thermal equilibrium between the system 



and the environment |||, or correlated initial states |T8l [19 



3 Stochastic Forces from the Influence Functional 

In this paper we will be interested in models in which the action (2.1) has a Hamiltonian of 
the general form 

H(a, q) = H(a) + H e (q) + Acr(a, a)e(q n , q n ) (3.1) 

n 

where A is a coupling constant and a and e are arbitrary functions of the system and envi- 
ronment variables. The simplification made in (3.1) is that system environment interaction 
is separable. This ensures that the effect of the environment on the system can be described 
by a single stochastic source. 
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Let us introduce the sum and difference system variables as 

E = - (cr(a, a) + <r(d, a')j , A = <r(a, a) — a(a', a'), 
and define the real quantities 



C n (ti, ...,t n ; E^, S tnjt J 



5 n W[E(s),A(s)] 



hJ 5A(ti)...5A(t n ) 



A=0 



(3.2) 



(3.3) 



where W = InjF. The notation of Ci(£i; E tl;t .] means Ci is a function of ti and also a 
functional of E between endpoints t\ and t«. Writing the influence action as a functional 
Taylor series and generalizing the notation to n variables we find that formally 



W[E( S ),A( S )] 



- ^2 d*i J*' dt 2 A(t 1 )A(t 2 )C 2 (t 1 , t 2 ; E tl )ti , E t2;ti ] 

+ ... + ^ (£) ^ dt 1 ...dt n A(t 1 )...A(t n )C n (t 1 , t n ; S tl;ti , .., S Wi 



+ 



(3.4) 



This form of the influence functional will turn out to be useful for its physical interpretation. 
From (2.9) and the propagator U given by 



U[a t ,u} = Y[Texp ~J ds(H e (q, s) + Acr(s)e(g n , <?„)) 



(3.5) 



it is observed that C n is of order A™ in the coupling constant. 

We can interpret the C n in (3.4) as cumulants of a stochastic force. Consider the action 



S[a{s)\ 



ds(^L(a, a, s) + a(a, d)£ (s 



(3.6) 



where £(s) is some forcing term. We want to consider the case where £(s) is a stochastic force 
with a normalised probability density functional V[cr(a, a); £(s)]. The probability density 
functional is taken to be conditional on the system history a(a, a). The action (3.6) generates 
the influence functional 



h 



.F[E,A] = ^exp 

= j Z^[£, E] exp £( S ) A( S )d S 



(3.7) 



where E and A are defined in (3.2). The essential point is that the influence functional 
(3.7) has the exact same form as the characteristic function of the stochastic process £. 
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Therefore given any influence functional A], we can interpret the C n , given by (3.3), as 
the cumulants of an effective stochastic force £(s) coupled to the system in a way described 
by the action (3.6). The probability density functional V[£, S] of £(s) can be obtained from 
a given influence functional by inverting the functional fourier transform (3.7). 

If we ignore all cumulants after the second order (the cumulants are of order A n ) we are 
making a Gaussian approximation to the noise. Although A is usually assumed to be small 
for the series (3.4) to converge, the formal expansion in orders of A is acceptable even if 
A = 1, as long as the deviations from Gaussian are small. With the Gaussian approximation 
we can write the influence functional as 



= jD£V[Z,T]exp^S IF [a,a',S] 



exp 



-S IF [a, a',f] 



(3.8) 



where 



P[£,E] = P exp(-| dt, J dt 2 i{t 1 )C^\t 1 M^t 1 ,t i ^t 2 ,tMt2)) (3.9) 



is the normalised functional distribution of £(s) and 



(3.10) 



We can use this effective action to obtain our semiclassical equation of motion which is 
given by 

5(S eff [a,a',t}) 



5A a (t) 

where A a = a — a'. We find it becomes 



A a =0 



(3.11) 



dL 



d da 



da dt da \ da dt da J 



I) (c&vtd+taj) - %(ci(t,*t, ti ]+m) = o (3.i2) 



where L(a,d) is the system Lagrangian and £(£) is a zero- mean Gaussian stochastic force 
with a correlator given by 



(3.13) 



Clearly both the noise and driving term are still dependent on the system history in a 
complex way in general. However we can further simplify things by expanding around a 
background a = a^. In this case we approximate the first cumulant by 



Ci(t;<7 titi ] = C^t-at^W^ + f dt'a(t')ti(t,t') + ... 

Jti 



(3.14) 
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d(t; (T tlti ] = d(t; a t , ti ]\ a=ab + / dt'a(t')fi(t,t') + ... 
where a = o — a b and 



n(t,t 



SCi(t; a tik ) 



fi(t, t 



5a(t') 
a _ SCi(t; <T t ,ti\ 



5a(t>) 



(7=a b 



u=a b 



(3.15) 



(3.16) 



(3.17) 



where we have assumed in (3.15) that fi(t,t') is antisymmetric as will be the case for our 
examples. The noise £(t) now has the correlator 



{t{tW))=C 2 {t,t'-cj tM ,cj t , M \\ (T -_ 



(3.18) 



These approximations greatly simplify (3.12). Our task is then to solve for the fluctuations 
a = a — af, subject to the initial condition a(tj) = a(£j) = 0. 



4 Influence Functional for Cosmological Backreaction 



In this section, following the methods of ||21|| , we will derive the form of the influence func- 
tional in terms of the Bogolubov coefficients in the transformation between the creation and 
annihilation operators of field amplitudes at different times. First we show how the dynam- 
ics of a general real scalar field in an expanding FRW universe can be described by a sum 
over quadratic time dependent Hamiltonians. Then we discuss the Bogolubov coefficients in 
terms of the squeeze parameters [31]. It also applies to the case of gravity wave perturba- 
tions whose two polarizations obey wave equations of the same form as a massless, minimally 
coupled scalar field (see |32| for details). 

The action for a free scalar field in an arbitrary space-time can be written as the sum of 
gravitation action S g and matter action S m of the form 



S 9 



l 2 p J d 4 x^g(R - 2 A) - 2l 2 p J d 3 xV^hK (4.1) 

S m = | J d*xV=9 (<T V, $ - (™ 2 + 2 ) + J d 3 xV^hK<5> 2 . (4.2) 

where I 2 = 1/(16ttG) and ^ = (n — 2)/4(n — 1) which in four dimensions d — 4 is equal 
to for minimal coupling and 1/6 for conformal coupling. Adding a surface term in the 
action proportional to K, the trace of the extrinsic curvature, is necessary for a consistent 
variational theory |33j in the treatment of the backreaction problem. 

In the spatially- flat Friedmann-Robertson- Walker (FRW) universe with metric 

ds 2 = a 2 (r ] )(d V 2 -Y^dx 2 ) (4.3) 
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R = 6a/ a 3 , K = 3d/ a 2 (where a dot denotes a derivative with respect to conformal time 
77 = / dt/a) we have 

S 9 = -Vf p J dri (6d 2 + 2Aa 4 ) (4.4) 



<i 4 x 



(X) 2 - £(x,i) 2 - 2(1 - QO-XX ~ rn 2 a 2 + (6£ - 1 



(4.5) 



Here x — is the rescaled field variable and V is the volume of the universe. From now 
on we will absorb l p by rescaling x an d a. 

We can expand the scalar field in a box of co-moving volume V (fixed coordinate volume) 



,. cos k ■ x + q r sin • x] 



(4.6) 



which leads to the Lagrangian 



2(1 - 6&) r g£g£ - ( k 2 + m 2 a 2 + (6& - 1 



(4.7) 



where k = \k\ and S m (rj) = J L(rj)drj. The canonical momentum is 
Defining the canonical Hamiltonian the usual way we find 



PS + (1 - 6^) T (« + + ( A: 2 + mV + 6&(6& - 1 



,<r2 



sk^k ^k^sk 



(4.9) 

where the sum is over positive k only since we have an expansion over standing rather than 
travelling waves. The system is quantized by promoting (p|, qZ) to operators obeying the 
usual harmonic oscillator commutation relations. In this way the dynamics of the field is 
reduced to the dynamics of time-dependent harmonic oscillators. (See [^l[ for details.) 

In the case of a free quantized scalar field coupled to a spatially- flat FRW universe with 
scale factor a(s) the action thus belongs to the general form 



S[a, q] 



ds 



L(a, a, s 



+ Y,{l m k( a ia)(Vk + h{a,d)q k q n -ujl{a,a)qfy} 



(4.10) 



By tracing out the scalar field we can obtain an influence functional and from which de- 
rive an equation of motion for the scale factor in the semiclassical regime. Here since we 
work explicitly in the semiclassical regime, the environment is quantum and gravity enters 
classically through the scale factor a. 
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We want to calculate the influence functional for this model. From (2.9) we see that it 
is formally given by 



where Ut is the quantum propagator for the bath mode in (4.10) with a(s) treated as an 
arbitrary classical time dependent function. We have derived this propagator before |H| and 
will only quote the results here. The result for a particular mode is (we drop the mode label) 

U[a tM } = S(r,<p)R(6) (4.12) 

where 

R{9) = e~ i6§ , S(r, 0) = exp[r(ie- 2 ^ - iV^)] (4.13) 

and 

A=j, ^ = ^-, B = tfa + 1/2. (4.14) 

S and R are called squeeze and rotation operators respectively. The parameters r, 0, 9 are 
determined from the equations 

a = —ig* f3 — iha (4-15) 

P = ihp + iga (4.16) 



where 
and 



a 



e- w coshr, (3 = -e^ 2 ^ sinhr (4.17) 



^-K^ 5 ^^)- (4 - l9) 

c is an arbitrary positive real constant that is usually chosen so that g = at tj. We must 
have a = 1 and /3 = at t,i so that the initial condition for the propagator is satisfied. The 
time dependence on g and h comes directly from a in (4.10). 

Applying (4.12) to (4.11) we find that the influence functional for a mode in an initial 
vacuum state is given by 

J^ k [a,a'} =Y,(n\S(r,<P)R(9)\0)(0\Ri(9')Si(r>, <f>')\n) (4.20) 

n 

where \n) are the usual number states. Using 

R(9)\0) = e- ld/2 \0) (4.21) 

we find 

F k [a,a'\ = E[(^|5(r,0)|O)(O|^(r',0')l^)]e- l(e - e ' )/2 . (4.22) 
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With 



S(r,<f>)\0) = (coshr)- 1/2 ^ 

n=0 



<(2nV 

-e^tanhr)" V - , \2n) 



2 n n\ 



and making use of 



we can show that 



E 



(2n)! n 

x 



[n 



\\2- 



1 



F[a, a'] = Yl 



1 



(4.23) 



(4.24) 



(4.25) 



* ^/«4a']afc[a] -/3 fc [a']/3|[a] 

This shows yet another way of deriving the influence functional in terms of the Bogolubov 
coeffients, in addition to the derivations given in |22|. 



5 Einstein-Langevin Equation 

From the Hamiltonian (4.9) we see that the system- environment interaction is separable for 
two cases: the massive conformally coupled field (for which a = a 2 in (3.1)) and the massless 
minimally coupled field (a = a /a) which also describes gravitons. For these two cases the 
results from Sec. 3 apply: (3.12) is the appropriate equation describing backreaction of the 
quantum scalar field on the metric. To derive the Einstein-Langevin equation we need to 
compute the first two cumulants given by (3.3) using the influence functional (4.25). 
The solution of (4.15) and (4.16) can be written as 

U[a t>t J = Texp (^—ij t ds u(s)j (5.1) 

■w ■ ( -m ) ^ 

The key to calculating the functional derivative of (5.1) is recognizing that we can always 
write U[a ti j.] = U[a tiT ]U[a T)t J. We therefore find 

^"W^^-W- (5 - 4) 

Making use of the formal expression for the time ordered representation of (5.1) it is easy to 
see that 

5\J\a tT ] r . r f , 5u(s) 



where 



and 
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5A(t) 



—i 



Su{s) \ 
as — , ; | U[a T)t J. 



5A(r) 



Substituting (5.5) and (5.6) into (5.4) we find that 

*U[a titi ] 



(5.6) 



(5.7) 



5.1 Massive conformally coupled field 

For the massive conformally coupled case we have a = a 2 and 
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(k 2 + m 2 a 2 )l 2 p c 



I 2 



, h 



(k 2 + m 2 a 2 )l 2 



+ 



(5.8) 



in (4.18-19). From (3.3) and (4.25) (we have reinstated the Planck length) we find the first 
cumulant of the stochastic force is 



l 2 m 2 +~ l 2 m 2 h + 



(5.9) 



k u k 

where q 2 = IJ^ [a 71>ri ^\q 2 tj[a r] ^ i } and the average is with respect to the vacuum. The propagator 

U is given by (4.12) with the Bogolubov coefficients determined via (4.15-16) with g, h given 
by (5.8). We will use this notation below as well. Similarly for the second cumulant we find 



74 4 



16 - 



ee4[(^+«,) 2 K'+^o' 



+ (0* + a;) 2 (a v , + fy) 
Applying (3.16) to (5.9) we find the dissipation kernel to be 



(5.10) 



MW) = -7]-EE[(tf>-(9}?} 

^ k 



ilphm 4 



E E i [(A, + «,) 2 K' + /^) 2 - (/?;; + «;) 2 (^ + /3,0 s 



8 "T c 



(5.11) 



Again we see the close relation between the noise and dissipation kernels. 

Using (4.4) and a = a 2 we find that the equation of motion (3.12) with the background 
approximation becomes 

a{rj) 



2 t o a(rj) 
a - -Aa 3 + 

3 6V1 2 



6VII 



(5.12) 



r f 7 ) i 

Ci(v,a 2 r,, Vi \\ai=al + / <Va 2 (?/V(?7, ?/) : 

where £ is a zero mean gaussian stochastic force with the correlator (5.10) evaluated on the 
background a^. 
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5.2 Massless minimally coupled case 

For the massless minimally coupled case, a — a/a, 



g = -i- + 



1 (Ilk 2 c 



a ' 2 \ c 11 



h 



we get 



2 I c II 



%h 



(5.13) 



dfa; (d/a)^] = ~EE((P9 + 9P).> = "tEEW " ( 5 - 14 ) 



where p is the canonical momentum from the Lagrangian (4.7) with m = £ = 0. For this 
so we see from (3.12) we must find C±. Taking the derivative of (5.14) 



case "o - 



da d do 



da dr\ da 



and using (4.15-16) and (5.13) (with c = Ilk) we find 



61(77; (a/a)^ Vi ] = h^2J2k[a*f3 v + a v (3*}. 



(5.15) 



For the second cumulant we find 

C 2 (ri,ri';(a/a) rm ,(a/a)rf, Vi ] = ^EE[((M + ?PWP9 + 9P)r/} + {{PQ + Qp) v '{pq + qp) v ) 

= t EE[K-^ 2 )«-^ 2 



+ (a* 2 -/3* 2 )(a 2 , -/3 2 ,) 
From (3.17) and (5.15) the dissipation kernel is given by 

m,i/) = -a EE *[(/?+ a ?) « 2 - fl?) + + O (4 - $ : 



a/a=(a/a)i, 



(5.16) 



(5.17) 



The equation of motion (3.12) with the background approximation becomes 



„_2 3 1 
° 3 A ° 12W|afa) 



61(77, {a/a) Vym )\ a/a={a/a)b + f drf ?i/i(?7,?/ 



12W2 a ( 7? ) 



. (5.18) 



We need to know the stochastic properties of £(77) given that £(77) is a zero mean gaussian 
stochastic force with the correlator (5.16) evaluated on a background. We can deduce this 
by integrating by parts the noise term in the effective action (3.4). We find that (relaxing 
the notation for C2) 



/ c?77i ^772 A (771) A (772)6*2 (771,772] = surface term 
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+ 



+ 



+ 



'If 

Vi 

dC 2 
din 



dC? 
drji 



dC 2 



dC2 

drji 



(m,Vf) T (vf) 



dr]i 



fVf 

dt x / dr] 2 T(r]i)T(7] 2 

•Irii 



d 2 C 2 (vuV2] 



(5.19) 



dr}xdr} 2 

where Y[j]) = J dt] A (77) = In a — In a'. The surface term will not contribute to the equation 
of motion but the last term of (5.19) shows clearly that £(t) corresponds to a zero mean 
gaussian stochastic force with the correlator 

The meaning of the middle term of (5.19) is more difficult to interpret. It vanishes only 
when the noise is stationary since we then have C 2 (r],r]'] = C 2 (i] — rf\. We will not discuss 
this term further since it will vanish in the example we consider next. Clearly though, its 
meaning will need to be considered for a study about nonstationary backgrounds. 



5.3 Backreaction of graviton fluctuations about flat space 

A simple case to study is a massless minimally- coupled field around a flat background 
(a = a). In this case ot{rj) = e~ tkr> and f3{rj) = 0. We see that in this case the first cumulant 
(5.14) vanishes. This should be compared to the massive field where the first cumulant is 
divergent around a flat background. The noise kernel (5.16) becomes 

fi 2 V r°° fr 2 V 
C 2 ( V -ri'} = — 2 dk k 2 cos[k( V - rj')] = ~S H {r, - rf) (5.21) 
Sin Jo oZtt 

where a prime on a function denotes a derivative taken with respect to its argument. From 

(5.20) we have 

C 2 ^ - rf] = ^ I dk fc 4 cos[fc(r/ - rf)) = ^""(V ~ v')- (5-22) 
The dissipation kernel (5.17) becomes 

fiy poo 

^ r] - ri ') = - dk k 3 coslkir] - rj')}. (5.23) 

l07T z JO 

The Einstein-Langevin equation (5.18) becomes 

& -l Aa3 - uvW) jC '">' Wf (n ~ r/) = (,,) - (5 ' 24) 

where £ is a zero- mean Gaussian force with the correlator (5.22). The solution of the 
Einstein-Langevin equations discussed here are beyond the scope of the present paper. We 
plan to consider these solutions in the future in the context of a general study into the 
dynamics of second order Langevin equations with non-local dissipation and colored noise 
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6 Summary 



Together with two related work [22|, P3|| , this paper seeks to establish a new framework for 
the study of semiclassical gravity theory based on the Einstein-Langevin equation. In 
the noise and fluctuation terms are identified from the closed-time-path formalism and the 
Einstein-Langevin equation derived for perturbances off the Robertson- Walker spacetime. 



In [2^] the influence functional method is used to derive an equation of motion for the 
anisotropy matrix of the Bianchi Type-I universe. Dissipation of anisotropy from particle 
creation in a quantum scalar field is seen to be driven by an additional stochastic source 
(noise) term related to the fluctuations of particle creation and shown to be a manifestation 
of a fluctuation-dissipation relation. In this paper, we have derived the following results: 

• By carrying out a functional Taylor series expansion on the influence functional we show 
how the successive orders measure the higher cumulants of noise in its most general 
(colored and multiplicative) forms, the lowest order truncation yielding a Gaussian 
noise. The second cumulant gives the autocorrelation function for the stochastic force 
(noise), which drives the Einstein-Langevin equation. 

• Using a general form for the Hamiltonian of a quantum field whose normal modes 
are coupled to a curved spacetime parametrically, we showed a new way to derive 
the influence functional in terms of the Bogolubov coefficients between the second- 
quantized operators of Fock spaces at two different times. This relation connects 
our new influence functional / effective action method with the traditional canonical 
quantization approach and thus incorporates the established body of knowledge in 
quantum field theory in curved spacetimes. 

• With the previous two results we were able to express the noise and dissipation kernels 
in terms of the Bogolubov coefficients. This connection offers a more transparent 
interpretation of the physical meaning of the many statistical mechanical processes 
such as decoherence and dissipation in terms of particle creation and related quantum 
effects. 

• We have also derived the form of the Einstein-Langevin equations for some well-studied 
cases of scalar fields in Robertson- Walker and de Sitter spacetimes. They form the 
starting points of the next stage of work, which is the solution of these equations for 
the analysis of fluctuations, instability and phase transition. We hope to report on 
these problems in future communications. 
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